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■ Abstract 



We calculate the gravitational radiation emitted by an infinite cos- 
mic string with two oppositely moving wave-trains, in the small am- 



plitude approximation. After comparing our result to the previously 
^ ' studied cases we extend the results to a new regime where the wave- 

{3J[). lengths of the opposing wave-trains are very different. We show that 

in this case the amount of power radiated vanishes exponentially. This 
means that small excitations moving in only one direction may be very 
^ ' long lived, and so the size of the smallest scales in a string network 

^ , might be much smaller than what one would expect from gravitational 

back reaction. This result allows for a potential host of interesting cos- 
mological possibilities involving ultra-high energy cosmic rays, gamma 
ray bursts and gravitational wave bursts. 



1 Introduction 

Topological defects are a prediction of most particle physics models that in- 
volve symmetry breaking and are therefore quite generic. They are formed 
when the topology of the vacuum manifold of the low energy theory is non- 
trivial The typical size of regions which acquire different vacuum expec- 
tation values depends on the type of phase transition but it can never exceed 
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the horizon size, due to causahty. The specific defect formed in turn depends 
on the symmetry of the groups involved, in particular on the homotopy group 
of the vacuum manifold. If the vacuum manifold contains disconnected parts 
domain walls are formed, if it contains unshrinkable loops strings are formed 
and if it contains unshrinkable spheres monopoles are formed. More compli- 
cated hybrid defects may be formed if there is more than one phase transition. 
For a review see 

The most simple phase transition where cosmic strings are produced is 

U{1)^1. (1) 

The vacuum manifold is U{1), which contains unshrinkable loops. The phase 
transition, occurring at energy r]s, traps the magnetic flux of the gauge fields 
associated with the U{1) symmetry into strings of mass per unit length 
/i ~ r/^. The resulting system is a network of long strings and closed loops. 
Numerical simulations and analytic studies show that this network quickly 
evolves into a "scaling regime" (see [|| and references therein), where the 
energy density of the string network is a constant fraction of the radiation 
or matter density and the statistical properties of the system such as the 
correlation lengths of long strings and average sizes of loops scale with the 
cosmic time t. The question of the size of the small-scale structure of cosmic 
string networks, however, has so far resisted concordance. Simulations show 
that most loops and the wiggles on the long strings have the smallest possible 
size, the simulation resolution. However, the prevailing opinion on this issue 
is that the size of small-scale structure also scales with the cosmic time t and 
its value is given by the gravitational back-reaction scale TGut, where F is a 
number of order 100 and G is Newton's constant. This possibility was first 
pointed out in |^. 

Cosmic strings have been a popular choice in the literature because un- 
like monopoles and domain walls they do not cause cosmological disasters. 
Indeed, they originally were considered good candidates for structure for- 
mation. Recent cosmological data, in particular the cosmic microwave back- 
ground power spectrum, is not consistent with density fluctuations produced 
only by strings. A combination of density fluctuations produced by inflation 
and strings is possible, and strings with small values of /i would have little 
gravitational effect. Cosmic strings are also good candidates for a variety of 
interesting cosmological phenomena such as gamma ray bursts Q, gravita- 
tional wave bursts |^ and ultra high energy cosmic rays |^. However, some 
of the predictions of these models depend sensitively on the so far unresolved 
question of the size of the small-scale structure. 
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Sakellariadou computed the radiation from a cosmic string with hehcal 
standing waves, which one can consider to be composed of excitations of 
equal wavelength traveling in the two directions. Hindmarsh ||8| calculated 
the radiation from colliding wiggles of different wavelengths on a straight 
string, in a low-amplitude approximation, and found that the radiation rate 
from small wavelength wiggles approached a constant as the wavelength of 
the wiggles in the opposite direction became large. However, Garfinkle and 
Vachaspati ^ showed that a wiggle moving on a straight string does not 
radiate. This poses a puzzle because an excitation of sufficiently large wave- 
length would seem to appear straight to a tiny wiggle passing over it. 

With this work we hope to shed some light on these issues. The following 
section contains a review of string motion and a derivation of the energy 
emitted in the form of gravitational radiation per unit solid angle. Section 3 
starts with a derivation of the power emitted from infinite planar excitations 
on a string with a subsequent specialisation to sinusoidal wave-trains in the 
case where the amplitude is small. We show that in the similar wavelength 
case the result of Hindmarsh ||8| applies and how an argument on the limit 
on the size of the small-scale structure consistent with |^] can be obtained 
from this result. We then demonstrate that this argument cannot generally 
be made and in particular show that when the wavelengths are very different 
the radiation vanishes exponentially. In section 4, we discuss these results 
and their relevance to the question of the small-scale structure in cosmic 
string networks and conclude with some remarks on the effect of these results 
on cosmological phenomena. 



2 Preliminaries 

When the typical length scale of a cosmic string is much larger than its 
thickness, 6s ~ rjj^, and long-range interactions between different string 
segments can be neglected, the string can be accurately modeled by a one 
dimensional object. Such an object sweeps out a two dimensional surface in 
space-time referred to as the string world-sheet. 

The infinitesimal line element in Minkowski space-time with metric ri^i, = 
diag(l, —1, —1, —1) is 

ds^ = r^^.dx^'dx' = r]f,,x%x';,dCd(:\ (2) 

where a = 0, 1 labels the two internal parameters of the string world-sheet 
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and x'^ = dx^/dC,"". One can then write the induced metric on the world- 
sheet of the string as 

(3) 



For an infinitely thin string we can use the Nambu-Goto action. It is pro- 
portional to the invariant area swept by the string, 



S = J dA = -n J d'^Cy 



-1, 



(4) 



where 7 = det(7ab) and /i is the mass per unit length of the string. The 
parametrisation-invariant energy-momentum tensor for a cosmic string is 
given by 

T'^^ix) = ^ / dW^r'^x^.S^^Hx - X(C)). (5) 



In the light-cone gauge the metric takes the form 



Jab 



1 

1 



ab 



1 



7(10, 



(6) 



because the arc element is ds'^ = 2x^u • x^ydudv with = x^{u) + x^{y), 
and x'^(n) and x'^v) are both null vectors. In this gauge the stress energy 
tensor (^) can be put in the form 



T^^'ix) =fi f dn(it;(« + x'^^xDS^^^^x - x{C)). 



(7) 



Far from a source localized in space and time, the total energy radiated in 
the form of gravity waves in the direction of k is |10| 

^ = 2Gj^ du;u.'{T'^^*{k)T,Ak) - - \Tli{k)\ }. (8) 

So what we need is the Fourier transform of the stress energy tensor 

1 



T^^^^k) = ^ I d'^xT''^{x)e 



ik-x 



(9) 



which is given by 



T^'^k) = d^xe'^ '^ j dudvixf^^x"^^ + x''^^xi^^)5'^^\x - x(C)). (10) 
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Figure 1: Two oppositely moving wave-trains on an infinite straight cosmic 
string. We ignore the effects of the kinks at the edges of the wave-trains 
because later we will be taking the length of the wave-trains to oo. 



fll) 



Using x^' = K(n) + bf'{v)]/2 yields for ([T 

r^^(A;) = —fif dudvia'^b'" + a"'b'^')e'''<''+'''>/^. 
Svr J 

Following this expression can in turn can be written as 



T^"'{k) = ^^l{A^'{k)B^k) + A^{k)B^'{k)) (12) 
Svr 



where 



-00 J —00 

We can then re-write the total energy in the direction of k (^) as 

'^-^0,^r^^H\AflBUlA'.Bf-lA.Bfy. ,14, 

3 Gravitational Radiation from Cosmic Strings 
3.1 Planar Excitations 

We consider planar string excitations moving on the z-axis in the two direc- 
tions as shown in Figure [T]. 

We are free to take 

a'''(n) = {l,f'{u),0,-^l-n, a^^iu) = (^n,/(n),0,- J ^1 - f'^du^ 

(15) 

where the sign of the square root is chosen such that u decreases in the 
positive z direction, which yields from (p^) 



A''{k)= r dua'^'{u)e^p- uju-k^f + k, f - f'^du 
J— 00 2 J 



(16) 
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exactly. The situation for B^{k) is completely analogous. 

Equation (H ) gives the energy radiated from a localized interaction. In 
Appendix we extend this to the power per unit length from an infinite 
periodic wave train. We find that the radiation is emitted in a discrete set of 
cones and at a discrete set of frequencies. Specifically, let be the average 

of Ki , 

^- = \- r duJl-f>^ (17) 
where is the wavelength of the wiggles of a^{u), and let 

t^a = ^ (18) 
Aa 

and similarly for and Kb. Then for each set of positive integers n, m 
satisfying 

(1 - A,)/(l + Ab) < \nKa/mK^,\ < (1 + A,)/(l - Af,). (19) 

we have radiation with polar angle 9 given by 

UKa - rriKb .^^x 
cosB = — , (20) 

riKaAb + rUKbAa 

at frequency 

" = ' A. + Ab ' (21) 
(so that uj + Aakz = 2nKa and uj — Abkz = 2mKb) with power per unit length 
per unit azimuthal angle 

dP 



dzd(p Aa + A; 



J2inKa + mKb)\An\'^\Bm\'^ (22) 



where 



'■n 



Aa J-\a,/2 

— / dua^ exp - 
Aa J-\ali 2 



<jJU + 



r ^ji-f'^du' - k^fiu) 



(23) 



If we fix the shape of a^, but vary the scale so that Aq becomes large, then 
Aa is fixed and — > 0. If we follow a single mode n, m as we go toward 
this limit, 9 reaches vr at a finite value of Ka, given by 

C = -^6(1-Aa)/(1 + Ab). (24) 
n 

For larger values of Ka, the given mode does not exist. 
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3.2 Sinusoidal waves 



3.2.1 Small amplitudes 



We will consider excitations whose amplitude is small compared to their 
wavelength, so that /' <C 1. This enables us to replace y/l — by its 
average, A^, in ([23|), to get 

1 /■^«/2 

~ ^ / dua'^'{u)expi[nKaU - k^f{u)/2] . (25) 

Aa J-\a/2 



In particular if we take the wave trains to be sinusoidal, namely, 

f{u) = eaCos{Kau), f{u) = — sin(K„ti) (26) 

with ea <C 1, and thus 

Aa^l-el/A, (27) 
for the time component we have 
1 r>^a/2 



1 f^a/-^ 

A^Ki — du expz [n^aU - Eak^ Sm{KaU)/2Ka] 

Aa J-Xa/'i 

which can be integrated to get 



(2^ 



Al - Jn ('^ ) (29) 



Ik a 

where J„ is the Bessel function of order n. The first component can be put 
in terms of Aq using integration by parts 

Ai{k) = (30) 

and the third component can be also be written in terms of as 

Al = Aa^o (31) 

so that 
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Figure 2: Discrete radiation cones given by (|2ClD and (|33|) for m = 1 and 
71 = 2,4 and 20 with = el = 0.1 and = Kb. Shading denotes the 
intensity of the radiation, with white being most intense. As n gets larger, 
the cone becomes narrow and the radiation gets concentrated in a narrow 
band around where the cone intersects the plane of the string. 



A similar analysis for B and using (|27D yields the expression for the power 
in terms of Bessel functions, 

dP Gfi^ ^. .(el in^Kl\ (el 



X 4 ('-^) Jii^-^] (33) 



2Kq 

As an example, Figure S shows some of the discrete radiation cones given by 



'el 


Ati?hI 








f ebkx \ 




\2kJ 



(P0| ) with power calculated from 



3.2.2 Similar wavelength v^^iggles and the back-reaction argument 

If Ka and Kb are similar in magnitude, and e^, <C 1, then the arguments of 
both Bessel functions will be much less than 1. In this case we can use the 
small-argument approximation for the Bessel functions [p]] , 

Jnix) - ^ (34) 
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and so we see that the dominant contribution to ( |33| ) comes from the n 
m = 1 mode, and the power can be written 

dP ^Gf^hlel 



dzdcj) 32 



'-(Ka + Kh). (35) 



Integrating over (p yields 



^ ^ ^^G'/^^^ae6('^a + l^b) (36) 

This agrees with the result of ||8|. The factor of 16 results from our choice of 
normalization in the definitions of A^^{k) and B'^(k). 

An argument that can be thought to follow from this result is the following. 
Consider short wavelength excitations in b interacting with long wavelength 
excitations in a. If we put 

Ha < Kb (37) 

in (^) we obtain 

-TT Gji e^e^— = -vr G^e„— (38) 



dz 8 ^ " " Afe 8 ^ Ab 

where 5^ ~ e^/x is the contribution of the short wavelength wiggles to the 
mass per unit length of the string. We expect the radiated energy to be 
taken from the small-scale wiggles, so we expect 

—t ~ ^^'^h ^''^ 

which has the solution 



with 



6fi oc e"*/^ (40) 
8Xb 



(41) 



If we assume that large wavelength wiggles exist with e ~ 1, we would 
conclude by the above argument that the minimum wavelength of wiggle 
that can survive until the present day is 

Amin ~ Gflto (42) 

where to is the present age of the universe. On smaller scales, wiggles are 
exponentially suppressed. 
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If we extend this analysis to a loop, the structure of the loop that enables 
it to be closed must involve features whose amplitude is comparable to their 
wavelength, so the condition e ~ 1 is always met. 

However, it turns out that ( ^6|) is never correct in the regime of (|37|) with 
ea ~ 1. Even 6^,66 <C 1 is not sufficient, because we might have, for exam- 



ple, kx/Ka 3> 1, and then the argument of Jn in (|33| ) would not be small. 
Explicitly, 

kx = u; sin(0) cos((/)) 



= A '^"^ A \/ {^bnKa + AaTUKhf - {uKa - niKhf (43) 

so kx is at most of order y/nniKaUb and thus we are in the small-argument 
regime whenever 

eln/Ka <. 1, elKa/Kb < 1- (44) 

Larger values of n and m could in principle make the arguments of the Bessel 
functions larger, but it turns out that in this situation the contribution is 
dominated by that of the lowest mode. 

One can write (|4^) in terms of the wavelengths of the sinusoids, Xa^b = 
'^'^/i'^a,b and their amplitudes, = €a^b/'<-a,b, as 



Aa,Ab^VXAb (45) 

which is invariant under Lorenz boosts along the string. Rather than merely 
requiring that the amplitude of each wave train be small compared to its own 
wavelength to have the simple case, we must require that each amplitude 
be small as compared to the geometric mean wavelength. Note that these 
conditions cannot hold if Xa 3> Xb and ea ~ 1- As an example of this, in the 
case of small "chiral" wiggles (those moving in one direction only) on a loop, 
the conditions ( ^ ) don't hold and we expect wiggles to survive on the loop 
longer than (p|). 

It should be noted that if both of the wiggles have about the same wave- 
lengths then an argument analogous to the one above can be made iQ. In 
this case wiggles in both directions contribute to the radiation power, and 
the extra energy in the string declines as 1/t. 

3.2.3 Chiral excitations 

In the following we take the small argument limit for the B part of 
namely e^Ka/^-b 1, or Ab <C \/ AaAf,, but do not use that approximation 
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for the A part of the expression. This corresponds to the situation where the 
opposing excitations have very different wavelengths. Thus ( ^3|) becomes 



For m = 1, the bounds (|19D become 



(1 - ^^)^>^ < ^ < ji + KW 



which using ( |27| ) can be written as 



X 8x 

- < n < (48) 

where 

X = (49) 

Let us first consider what happens to the n = 1 mode as Ka approaches the 
limit in (p^). In this limit, 6 becomes vr, while uj and approach fixed, 
finite values. Thus kx/na — > 0, we can take the small argument limit of the 
Bessel function and the power in this mode approaches a finite value given by 
(|36|) . When Ka exceeds the limit of (p4|), the n = 1 mode vanishes suddenly, 
causing a discontinuity in the emitted power as a function of Ka- Similarly, 
if we let ea approach the limit given by ~ 8Ka/Kh, we see that the power 
in the n = 1 mode is again given by (^), which increases with until it 
suddenly vanishes. 

As Ka decreases, the modes successively vanish (although only the energy in 
the n = 1 mode does so discontinuously) , but the remaining modes that are 
far from the limiting wave number grow larger. To see the overall effect, we 
consider the case where e^Kh/Ka ^» 1, or in terms of the amplitude Aa ^ 
VAaAfe, so that many different modes contribute to the total radiation. 

In Appendix |b| we use the large-order approximation for Bessel functions 
to compute the power ( ^ ) in this regime. We find that the power declines 
exponentially with x, 

— ~ Gi?elKae-^^ (50) 

where a « 0.07867. Thus there is essentially no radiation from wiggles of 
significant amplitude and very different wavelengths. 
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To compare the analytic approximation ( [50| ) with our original expression 
(^6|) , in Figure § we show a plot of ^ versus x = e'^Ki,/Ka with = .1, 
eb = .05 and Kq fixed. The dashed line is 

with a best-fit value /3 = —5.69, and a = —0.07867 as above. The points 
are given by the sum over n and numerical integration over (p of (|4^). In the 
large x regime, we see that the agreement with the analytic approximation 
is excellent. 

The exponential suppression we have just found is valid when the Lorentz 
invariant conditions 

Ab < ^KAb, Aa > (52) 

are satisfied. It is interesting to consider the effect of a Lorentz boost that 
makes both wiggle wavelengths similar. This boost makes the condition (|5^) 
for Aa read 1- In this case we could not have performed the calculation 
at all because this is not the regime in which (^) is a valid approximation. 
It must be true, however, because of Lorentz invariance, that the radiation 
is also exponentially suppressed. Although at first glance this may seem 
puzzling, it is not too difficult to understand this result: In this limit, even 
though the wavelengths are similar, the amplitude of one of the wave trains is 
much larger than the other {Aa ^ Ab) so that the larger wave-train appears 
to be traveling on an almost straight string and would not be expected to 
radiate 



4 Discussion and Conclusions 

We have calculated the power radiated in the form of gravitational waves 
from oppositely moving excitation wave-trains in the small amplitude regime 
when the wavelengths are comparable and when they are very different. 
We have shown that small opposing excitations on strings will only radiate 
significantly if their wavelengths are comparable. In this case we can always 
use (|^) for the power and the extra energy in the form of wiggles decreases 



as 1/t as discussed in the literature |12| . 



If, on the other hand, the wavelengths are very different then the grav- 
itational radiation can be suppressed. If the amplitude Aa of the long- 
wavelength excitations is greater than the geometric mean wavelength, \/Xa\, 
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Figure 3: Plot of In (^ q^a,. vs. x = e^Kb/Ka comparing our analytic 
approximation ( ^^ ) with f3 = —5.69 and a = —0.07867 (dashed line) with 
(H6|) calculated numerically for = .1, £{, = .05 (points). 



then the radiation rate is exponentially small. Thus unless the amplitude of 
the long-wavelength excitations is very small compared to their wavelength, 
there is little interaction between modes of very different wavelengths. 

The situation where the two wavelengths are very different occurs in loops 
with chiral wiggles, i.e., when the wiggles are only on the right- or left- 
moving excitations. In this case because the radiation vanishes exponentially 
the wiggles may be very long-lived. We therefore expect that if the loop is 
originally in a stable (non self-intersecting) trajectory it will radiate and 
shrink until it self-intersects, at the latest when the string loop is the same 
size as the wiggles living on it. This may take a long time because the tension 
on the strings is reduced by the presence of wiggles and the motion may be 
slow [p^ ]. 

When an excitation travels on an infinite string, it will eventually be exposed 
to every possible sort of excitation going the other way, and will thus eventu- 
ally lose its energy by gravitational radiation. However, when an excitation 
travels on a loop, it will encounter the same oppositely moving excitations 
over and over again. It is possible that when a loop is formed there will be, 
by chance, an excess of (say) right-going energy over left-going energy over 
a broad range of wavelengths. In this case, energy will be radiated until the 
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left-going excitations have been eliminated and some energy still remains in 
the right-going ones. Thus the string may become chiral with respect to a 
certain range of wavelengths, and the chiral wiggles in this range may then 
become long-lived. This process is quite complicated, and we will not try 
to analyze it here. We merely note that certain initial conditions may lead 
to long-lived chiral wiggles over certain wavelength ranges, and we briefly 
investigate the consequences of such excitations. 

Almost all proposals for cosmic-string-induced astrophysical phenomena are 
related to cusps. These are regions of the string that achieve extremely 
high Lorentz boosts Q. Therefore here we will focus on the effects that the 
existence of small wiggles may have on cusps. 

Near a wiggly cusp the effective mass per unit length is 

I2,s = + (53) 

where e is the amplitude to wavelength ratio of the wiggles and 7 is the 
(transverse) Lorentz boost. If we think of the string as having a thickness due 
to the wiggles then before the cusp can be formed the wiggles will typically 



overlap when the gamma factor reaches the value ||14| 



7o ~ ^JL/d (54) 

where L is the size of the loop on which the cusp appears, d ~ e/^ is the 
effective thickness due to the wiggles and k ~ 1/A is the wiggle wavenumber. 
The cusp may also not form at all if there is sufficient back-reaction from the 
wiggles to deviate from regular Nambu-Goto motion of the string, namely 



when the gamma factor in (53) reaches the value 

lb ~ (55) 

Which of the two processes is the dominant one depends on which one of the 
two gamma factors is the smallest, since it will be the one reached first by 
the cusp. 

It follows from this that there are two distinct wiggliness regimes. If 7;, > 70, 
or LKe < 1, there is not much change in the effective mass of the string and 
the wiggles will not significantly affect the string motion. However, before the 
cusp can be formed, because of the overlap of the wiggles, we will typically 
have a self-intersection that chops off a loop of size I ~ \/eL/K which in turn 
will fragment into about \J nLje loops of the wiggle size e/k. 
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Figure 4: Intersection on the unit sphere of a wiggly a' and a straight b'. 
The dotted hue resulting from averaging out the wiggles in a' makes an angle 
e with b'. 



If, on the other hand, 70 > 7?,, or Lne > 1, we are in the situation where 
the back-reaction is important. Much work has gone into understanding the 



behaviour of strings in this regime |15|. Typically one averages the string 



over scales larger than the wiggle size and the effect is to increase the effective 
mass density and decrease the tension of the string. The situation is then 
remarkably similar to that of a superconducting string with a chiral neutral 
current and we expect the cusp in this case to be smoothed out and self- 
intersections to occur near it some of the time as discussed in ||T6|. 



We can try to estimate the number of cusps formed in both cases. If we 
consider the section of the unit sphere where a wiggly a' and a straight b' 
intersect, as in Figure ^, we can see that the number of crossings must be 

A^c ~ e cot 61/1 Aa'l (56) 

where e cot 9 is the length available for crossings in the section where the a' 
and b' intersect on the unit sphere and |Aa'| is the distance between peaks 
on the unit sphere. Typically 

|Aa'| ~ |a"|Ao- ~ \/L (57) 

where A is the wavelength of the wiggles and L is the size of the loop. Taking 
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Figure 5: A "cusp" in a wiggly string with Lue > 1. The wiggles have led the 
overall shape of the string to be smoothed out. Instead of one large cusp, 
we expect a large number of small cusps. 

cot ~ 1 gives 

Nc ~ LKe (58) 

for the number of cusps. It is interesting to note that this is also the quantity 
that determines whether we are in the overlap or back-reaction regime. This 
means that in the overlap case there will still be one cusp, just as in the case 
where there are no wiggles. As we have argued, however, self-intersections 
due to overlap will take place before the cusp can be fully formed. In the 
back-reaction dominated case instead of a large cusp we will have about Lue 
small cusps on the smoothed out curve. In Figures ^ and ^ we plot pictures 
of Burden loops onto which we have superimposed chiral helical wiggles in 
the LKe > 1 case. 

From this analysis it appears that wiggly cosmic string cusps may solve some 
of the problems currently faced by cosmic-string-induced astrophysical phe- 
nomena. Wiggly cusps could produce substantially more ultra-high energy 
cosmic rays because of the presence of many small cusps. Gamma ray bursts 
may not be repeated because of self-intersections near the cusp. Gravita- 
tional wave bursts may be more frequent since, like cosmic rays, most of the 
energy comes from a small area around the cusp. Finally, wiggles may also 
lead to fine structure in gravitational wave and gamma ray bursts. 

Although some of these require different wiggliness regimes and would ap- 
pear mutually exclusive in fact they are not; we can "have it both ways" 
because as we have shown modes of different wavelengths barely interact 
and as a consequence chirality can be independently established at different 
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Figure 6: Another possibility for a "cusp" in a wiggly string with Lk€ > 1. 
In this case the smoothing effect of the wiggles has led to a self-intersection. 

wavelengths. 
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A Periodic sources 

In this appendix, we go from (p^, which gives the energy radiated from a 
localized interaction, to the power per unit length for an infinite periodic 
source. To do this, we consider a wave train of length L in a and another of 
length / in 6, and then let L and / go to infinity. The resulting energy will be 
proportional to LI, as shown below. We will ignore any contributions coming 
from the ends of the wave trains, because these will not be proportional to 
LI. 

Thus we will take / to be an odd periodic function in u, for Na periods of 
length Xa (with L = NaXa) centered around u = 0, and / = outside this 
range. Since / is odd, all components of will be real and (U) becomes 

dn 167r2 Jo I I I I y J 
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Because we are interested only in terms proportional to LI, we exclude the 
contribution of the straight parts of the string to the integral in ([T^), and 
write this equation as 

j=-(Ar„-l)/2'^(-'~^/2)'^'' 

Now a' is periodic, and so a{u + jXa) — ci{u) is a constant independent of u. 
Thus 

i=-(Afa-l)/2 -'-A<,/2 



where 



= UJ + k^Aa = io{l + AaCosO) (62) 
Aa = ^ d^x Jl - /'2. (63) 



with 



The sum in (|6lD is 

^"^"V^' c^i^ = sin(iVaa: J2) 

with Xa = XaK+/2. This function has peaks of height Na whenever Xa is 
an integer multiple of 27r and in the vicinity of these peaks the function is 
effectively NaSmc{NaXa/2Tr). However, because what appears in (|59| ) is the 
square of A'^{k) what we actually need to evaluate is sinc^ {NaX a /2tt). In 
the large Na limit Na s'mc^ [Nau) —>■ 6 {y) so we can write its square as 

sin2(iV,x/2) ^^^ g 5(^ - n) = 2.,iV, £ 5(K+ - 2n.,). (65) 
sm x/z zlzvr) 

' ' n=— oo ^ ' n=— oo 

A similar set of expressions can be written for B^^{k) which allow (^ to be 
re-written as 

^ -Ka^bA^aiVfe - 2nKa)5{K_ - 2mKb)A^B^ (66) 



dujdCl 47r2 



18 



where 

^=/" dna'V'="(")/2^ e= / ' dvb'''e'''-^^^^/\ (67) 

All the variables with the b subscript have the obvious meaning and = 
Lo{l — Aft COS 6*). The sum in ([66|) is performed for positive values of n and 
m only since and K- are nonnegative. 

We can now use the 6 functions to perform the integral in (|59|) over oj and 
to integrate over 6. We do this by a suitable change of variables, namely, 

io = (AfeK++A,ir_)/(A„+Afe), cose = {K+-K^)/{A,K++AaK^)) (68) 

so that 

(kodcos9 = dK+dK_/{AbK+ + AaK_) = dK+dK_/{uj{Aa + A^)). (69) 

Upon integration we get 

dE Gfi^ 



where the vector 



A , A LlJ2^MBm\' (70) 



Aa J-\a/2 

and there is a similar expression for B!^. Dividing (70) by the volume of the 
world sheet where the interaction occurs, yields the power per unit 

length per unit azimuthal angle (f) 

- ^'^ 5^(nKa + mKb)|AP|5mP (72) 



dzdcl) Aa + Afe 



Keeping the above in mind we can see that the significance of the appearance 
of the 6 functions is that for periodic sources the radiation is only emitted in 
discretized directions and frequencies. Explicitly, for a set of given Aq, A;,, 
Ka and Kb we have that 

cose= . (74) 

nKflAfe + rriKbAa 
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This discreteness in frequencies and angles arises from the temporal and 
spatial periodicity of the system. 



If we had = = 1, then the right hand side of { \14\ ) would always lie 
between -1 and 1. However, for waves of finite amplitude, this is not the case 
and thus not all modes are present. This limitation restricts the sum over n 
and m to satisfy 

(1 - A„)/(l + Ab) < \nKa/mKb\ < (1 + A,)/(l - A^). (75) 



B Large order approximation 

In this appendix, we compute an analytic approximation to the radiation 



power ( [4q ) using the Bessel function approximation for large order |11| 



gn(tanh a— a) 

Jninsecha) —^^^=^^= (76) 
V27rntanha 



which is valid when tanha > 1/n. 
Using (|4^) we can see that 



kl _ cos^((/)) 



{iA,n + Aax/elf-{n-x/elf} 



which allows us to express the square of the argument of the remaining Bessel 



function in (46) as 



k^el cos'^{<f>) ] ^ e\<i\ 2I 2 



and therefore 



Bearing in mind ( [7^ ) we take 

sech a = w = y/n k, cos{(l)){Ax/n - x^/2n^}^/^/4 (80) 



20 



so that we can write 



dP Gfi^ele?^^ , / 1 1\ , 



dzdcj) 8 2 

In our case it is easy to see that w has a maximum at hm = x/A and the 
value of the function at this peak is Wm = cos{(j))/y/2. This means that w 
is at most l/\/2 and thus tanha is at least l/V^, and we see that ( |76| ) is 
indeed a good approximation. 

We will now perform the sum over n. To do the calculation analytically we 
must first approximate the sum by an integral, which we can do provided the 
values of the function are sufficiently close to one another in the region in 
which we are performing the sum. Secondly, we will assume that the largest 
contribution to the Bessel function comes from a maximum and that we can 
approximate that integral by multiplying the height at the peak times the 
width of this peak at the maximum 



dP 

dzd(f) 



J2 Jninw) ~ J dnJl{nw) ~ J^ax^^. (82) 



This will be a reasonable approximation provided the function is sufficiently 
smooth around the maximum. Below we show that this condition is satisfied. 



We can take 

where 
If we let 



^'^^c'—il- — ] (83) 
2n V 8n ' ^ ' 



cos 



(l)/V2 (84) 



P=^-l (85) 
then 

w^ = c^{l-p^) (86) 

which obviously has the maximum value at p = 0. To find the maximum 
of the integrand of ( |82|) we use the very rough approximation 



In Jn{nw) ~ n(tanha — a) = n{\/l — w"^ — sech ^ w) . (87) 

Now if the factor of n were not present, we would reason that since w = sech a 
is a decreasing function of a, and tanha — a is also a decreasing function 
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of a, then the right hand side not counting the n is an increasing function 
of w, and (|87| ) will have its maximum where w has its maximum. Since the 
right hand side is negative, and also a decreasing function of n, we can infer 
that the actual value of n that maximises ([87|) is actually less than the one 
that maximises w. We can thus write 

" 4(1 +p) " MiTvT^^j^) ^^^^ 



and 



Setting the derivative of (|89|) to zero gives 

1 — — Vl — w"^ sech~^ w 1 — w'^/c^ + \/l — w^/c^ 
— w w/c^ 

This equation can be solved (although not analytically) for Wm, the value 



, . X X ^\ — w"^ — sech ^ w , , 

In J„(nu;) ~ j ^ ^ ,^ • (89) 

4 1 + ^1 — I& 



(90) 



of w that maximises (87) for a particular c and from it, using (|8q), one can 
obtain n^, the value of n at the maximum of ([89|). Note that both values 
will be in general different from n^, and the maximum of n and w at 
the maximum of the function w. Figure ^ shows a plot of Wm-, the maximum 
of (p^), as a function of c = cos((/>)/\/2- As expected the value of w where 
(|87|) is at a maximum is slightly less than the maximum value of w. 

In order to estimate the width at the peak we can again use (87). However 



in this case we can use the width at the maximum of the function i.e., 
at p = 0, or Um = x/A, because it is sufficiently close to the maximum of 



Jninw). This simplifies the calculation considerably. Prom (87) one can see 
that ^ 

— In J ~ — n tanh^ a (91) 
da 

and we can compute the width using the Taylor expansion 

In J ~ In 

Taking In J — In Jq ~ — 1 we find 



In J ~ In Jo ~ Aan tanh^ a. (92) 



Aa (93) 

n tanh a 

We need to express our answer in terms of An however so we take 

Aw ~ Aa = sech a tanh qAo = .. (94) 

da n^Vl - 
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and 

A- = ^4^, (95) 
because the first derivative is at the maximum, yielding 

/ ~ \ 1/2 

An=U=^ ^n]l\ (96) 



As advertised above, since fwn oc (e^K^/Ka)^/^ is large in the small Ka limit, 
the function around the maximum is sufficiently smooth that we can now 
perform the sum by approximating it as an integral and further approximat- 
ing tiiis integral by multiplying the value of the function at the peak by the 
width near the peak An. 

For the power per unit length per unit azimuthal angle this procedure yields 
where rim and Wm are given above. 

The integral over cf) can be approximated in an analogous way, taking the 
value of the Bessel function at the peak, which is at ^ = 0, and finding the 
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width at this peak. Evaluating the function at the maximum is straightfor- 
ward, one simply needs to solve ( pO|) for cos </> = 1, or equivalently c = l/\/2. 
This procedure yields wq ~ .691584 and from (B8), uq x/4. 83352. 

In this case the width can be found by using In J ~ f{wm) as given by ( |89| ) 
and we can see that 

dcj)"^ dwm dcp"^ 
because dwm/dcj) = 0. Using In J — In Jq ~ —1 again yields 

\2 Jf J2„ 



(98) 



{MY df d'Wm 

— X 

2 dwm dcfp' 



(99) 



and because / oc x, ignoring factors of 0(1) we have that 

A0~x"^/2. (100) 



These considerations yield for the power per unit length 

dP 

dz 



dP 

— ~ GfYelelinoKa + Kb) (mwo) AnA(j) 



Gfi^elel{noKa + k^) (mwo) (101) 



Expressing this equation instead using (|76|), with inverse hyperbolic functions 
written out explicitly, and ignoring factors of 0(1) yields 

which can be written more intelligibly as 

dP 

— ~ GfihlKae-^^ (103) 
where a ^ 0.07867 can be calculated from wq and uq. 
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